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ABSTRACT 

We calculate the 10 x 10 two-loop anomalous dimension matrix to order 0{aas) 
in the dimensional regularization scheme with anticommuting 75 (NDR) which is 
necessary for the extension of the AS" = 1 weak Hamiltonian involving electroweak 
penguins beyond the leading logarithmic approximation. We demonstrate, how a 
direct calculation of penguin diagrams involving 75 in closed fermion loops can be 
avoided thus allowing a consistent calculation of two-loop anomalous dimensions in 
the simplest renormalization scheme with anticommuting 75 in D dimensions. We 
give the necessary one-loop finite terms which allow to obtain the corresponding 
two-loop anomalous dimension matrix in the HV scheme with non-anticommuting 
75- 
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1 Introduction 



The inclusion of next-to-leading QCD corrections to the effective low energy Hamilto- 
nian for non-leptonic decays requires the calculation of the relevant two-loop anomalous 
dimension matrices. In the case of the AS = 1 Hamiltonian there are ten operators 
Qi,i = 1, . . . , 10 to be considered: current-current operators {i = 1,2), QCD penguin 
operators (i = 3, . . . , 6) and electroweak penguin operators (i = 7, . . . , 10). Consequently, 
one deals with 10 x 10 anomalous dimension matrices. 

Because of the presence of electroweak penguin operators with Wilson coefficients of 
0{a), a consistent analysis must involve anomalous dimensions resulting from both strong 
and electromagnetic interactions. Working to first order in a but to all orders in as, 
the following anomalous dimension matrix is needed for the leading and next-to-leading 
logarithmic approximation for the Wilson coefficient functions, 



The one-loop matrices ■y^^'^ and 7^°-' have been calculated long time ago []T|-|]TT|. The 
2x2 submatrix of the two-loop QCD matrix 7^^^^ involving current-current operators Qi 
and Q2 has been calculated in ref. [l^, |l^. Recently, we have generalized these two-loop 
calculations to the penguin operators Qi, i = 3, ... ,10 ||T^, |T5[, so that the full matrix ■y^^^ 
is also known. The purpose of the present paper is the calculation of 

The calculation of -^^l^ proceeds in analogy to 7s^^-' . In fact all the singularities calculated 
in ref. [|T3[ can be used here so that our main task was the calculation of the relevant colour 
and electric charge factors in two-loop diagrams involving one gluon and one photon. Yet 
as will be seen the structure of basic expressions and of the results differs from the one found 
in the pure QCD case, because now the electric charges of quarks matter and the flavour 
symmetries present in ref. [15] are broken. In fact this breakdown of flavour symmetry is 
the origin of the existence of electroweak penguin operators. 

In order to make the comparison with the calculation of 7g^^-* as easy as possible, we have 
organized the present paper in a similar way as it was done in ref. |1^. We present here 
however only the results in the naive dimensional regularization with anticommuting 75 
(NDR). The consistency of the NDR and 't Hooft-Veltman (HV) scheme for the calculation 
at hand has been demonstrated in ref. [|15[. There we have developed a method which allows 
to avoid a direct calculation of penguin diagrams involving 75 in closed fermion loops 
which are ambiguous in the case of the NDR scheme. Our method can be generalized to 
the present case so that an unambiguous calculation of 7^^ can be accomplished in the 
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NDR scheme. This is gratifying because this scheme is certainly the most convenient and 



contrary to the HV scheme [|1^, |T3|, satisfies all Ward-identities in the context of the 
minimal subtraction scheme. 

Our paper is organized as follows: In section 2, we recall the explicit expressions for 
the operators Qi, and we classify the one- and two-loop diagrams into current-current 
and penguin diagrams. We also discuss the basic formalism necessary for the calculation 
of In section 3, we recall the matrix ■yjP\ In section 4, the calculations and results for 
two-loop current-current diagrams are presented. In section 5, an analogous presentation 
is given for two-loop penguin diagrams. In section 6, we combine the results of the previous 
sections to obtain 7^'' in the NDR scheme. We discuss various properties of this matrix, 
in particular its Large-N limit. Section 7 contains a brief summary of our paper. In 
appendices A and B, explicit expressions for the elements of the 10 x 10 matrices 7^°-' and 
■y^l^ for arbitrary numbers of colours (A^) and fiavours (/) are given. In appendix C, the 
corresponding results for 7^^^ in the case of = 3 are presented. 

2 General Formalism 
2.1 Operators 

The ten operators considered in this paper are given as follows 



Qi = 


(Sa'U/3)v_A ('"/3'^a)v-A ' 


Q2 = 


(sm)v_a (^C?)v-A ' 


Q3 = 


(srf)v_A (^?)v-A ' 

q 


Q4 = 


(Sarf/3)v_A (^/35'a)v-A 
Q 


Q5 = 


(^^)v-A Y (^^)v+A > 
<? 


Qe = 


(Sa(i^)v_A Y ('ZS'?a)v+A 

q 


Qi = 


3 

^ (^C^)v-A J2 ^q i^l)v+A 

^ q 


Qs = 


3 

^ q 


Q9 = 


3 

9 {sd)y_p^Y.^q fe)v-A 

^ q 



(2.1) 



a;V+A ' 
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3 

where a, j3 denote colour indices = 1,...,A^) and are quark charges. We omit 

the colour indices for the colour singlet operators. {V ± A) refer to 7^(1 ± 75). This basis 
closes under QCD and QED renormalization. 

At one stage it will be useful to study other bases, in particular the basis in which the 
first two operators are replaced by their Fierz conjugates, 

Qi = (s(i)v-A (^m)v-a ' 

Q2 = (Sarf^)v_A (^/3^c^)v-A ' (2-2) 

with remaining operators unchanged. In fact the latter basis is the one used by Oilman and 
Wise 0. We prefer however to put Q2 in the colour singlet form as in eq. (|2.1|), because it 
is this form in which this operator enters the tree level Hamiltonian. Let us finally recall 
that the Fierz conjugates of the (V — A) ^ {V + A) operators Qi, i = 5, . . . , 8 are given by 

Qe = (^l^r) (qudi^) , 

q 

Qs = -12^e,(sLgR)(gRrfL) , (2.3) 

with similar expressions for Q5 and Qr- Here qn^i = |(1±75)q'- The Fierz conjugates of the 
{V — A)®{V — A) penguin operators Q3, Q4, Qg, and Qio to be denoted by Qi, i = 3,4, 9, 10 



are found in analogy to (2.2). 



2.2 Classification of Diagrams 

In order to calculate the anomalous dimension matrices 7^°-* and one has to insert 
the operators of eq. ( p.l|) in appropriate four-point functions and extract 1/e divergences. 
The precise relation between 1/e divergences in one- and two-loop diagrams and one- and 
two-loop anomalous dimension matrices will be given in the following subsection. Here, let 
us only recall that insertion of any of the operators of eq. ( p.l|) into the diagrams discussed 
below results into a linear combination of the operators Qi. The row in the anomalous 
dimension matrix corresponding to the inserted operator can then be obtained from the 
coefficients in the linear combination in question. 

There are three basic ways a given operator can be inserted in a four-point function. 
They are shown in fig. |1|, where the dot denotes the interaction described by a current in 
a given operator of eq. ( p.l|) and the wavy line denotes a gluon or a photon. 
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We will refer to the insertions of fig. |l](a) as "current-current" insertions. The insertions 
of fig. |I](b) and (c) will then be called "penguin insertions" of type 1 and type 2 respectively. 

The complete list of diagrams necessary for one- and two-loop calculations is given in 
figs. 0-^. At one-loop level one has three current-current diagrams (fig. 0) to be denoted 
as in ref. [|13| by Di - D^, and one penguin diagram (fig. ^ of each type to be denoted 
by Pq^^ and Pq^^ for type 1 and type 2 insertions respectively. At two-loop level there 
are 21 current-current diagrams shown in fig. to be denoted by D4 - D2A , and 12 
penguin diagrams of each type shown in fig. ^ to be denoted by P\^^ ^-PiP and P^"* 
respectively. The notation is similar to ref. [|l5l, but this time the wavy lines denote a 



gluon or photon so that complete 0{aas) is obtained. The omitted diagrams Pq and P7 
having triple boson vertices do not contribute here. The penguin diagrams which have 
no l/e divergences and some examples of penguin diagrams which vanish identically in 



dimensional regularization are given in figs. 6 and 7 of ref. |T^. It is needless to say that 
all possible permutations of gluons and photons have to be considered as well as left-right 
reflections. In the case of current-currrent diagrams also up-down reflections have to be 
considered. 

2.3 Basic Formulae for Anomalous Dimensions 

The anomalous dimensions of the operators Qi, calculated in the MS scheme, are obtained 
from the l/e divergences of the relevant one- and two-loop diagrams with Qi insertions 
and from the l/e divergences in the quark wave-function renormalization. Let us denote 
by Q a column vector composed of operators Qi. Then 

^{g,e) = Z-'fx-^Z, Qb = ZQ, (2.4) 

— * 

where Qb stands for bare operators. Working in D = 4 — 2 e dimensions, we can expand 
Z in inverse powers of e as follows 

00 1 

Z = i + Y.-Mg,e), (2.5) 

k=i ^ 

where g and e are the QCD and QED renormalized coupling constants. 
Inserting (|2.5|) into (^.4|), one derives a useful result 



2dZi{g,e) ^dZi{g,e) 



iM^-2,'^-2e'^. (2,6) 

Let us next denote by T^'^\Q) and T^^\Q^) the renormalized and the bare four-quark 
Green functions with operator Q insertions. Strictly speaking r^'^^{Q) and r^^\Q^) are 
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matrices, because the insertion of a single operator in a given diagram results in a linear 
combination of operators. 

At the one-loop level Fg'* is obtained by evaluating the diagrams of figs. 2 and 3. At 
the two-loop level it is found by evaluating the diagrams of figs. 4 and 5 and subtracting 
the corresponding two-loop counter terms. Next 

T^'\Q) = Z^Z-'r'i\Q^) , (2.7) 

where is a matrix which represents the renormalization of the four quark fields on 
external lines. In the pure QCD case this matrix was diagonal. The 0{a) and 0{aas) 
terms in this matrix are however non-diagonal, because in the case of penguin operators 
the dd, ss, and bb parts are renormalized differently from uu and cc parts. We next expand 
^ and Z^ in inverse powers of e as follows 



Z^ = i + J2-Z^,ki9,e), (2.8) 
k=i ^ 

T^^\Q,) = 1 + E -ZPig, e, Q.) + finite , (2.9) 



where 



k=i ^ 



^(0)^ e^g^ ^(1) 



e) = ^^^^ + j^,Z:;> + . . . , (2.10) 

10 / 2 2 2 \ 

zP{9,e,Q.) = E [j^,{d,h + J^M^h + •••], (2.11) 

where the dots stand for the pure QCD case which has been already considered in ref. |]15 
and higher order corrections. 

Demanding F'-^-' to be finite, we find Zi{g,e) and using ( 2.6 



(7f)^, = -2[(zf),, + (rf0.,] , (2.12) 
(7ieO„. = -^[izfh, + {d,),,] . (2.13) 



The matrices Z^p^ and Z^^ are given in section 2.5. 

2.4 Renormalization Scheme Dependence of 7(g) 

The two-loop anomalous dimension matrix depends on the renormalization scheme for 
operators and in particular on the treatment of 75 in D 7^ 4 dimensions. This is signaled by 
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the scheme dependence of the finite terms fg and fe in the renormahzed one-loop matrix 
elements 

(Q) 



as ^ a ^ 



where {Q^^^) denotes tree-level matrix elements. 

If Za and Zh are two renormalization factors of 
ization schemes a and 6, then 



1 + -^Afs + — Af, 



(QW), (2.14) 
corresponding to two renormal- 

(2.15) 



where 



Afs = (fs)b - ifs)a ■ 



Af, 



[Te b - [Te a ■ 



(2.16) 



Using eqs. ( |2.15|) and (p.4|) one finds the relation between calculated in schemes a 
and b. 

{lil^), = (7^1 + [Af.,7f)] + [Afe,7f)] . (2.17) 

This relation is very useful as it allows to test the compatibility of the two-loop ano- 
malous dimensions calculated in different renormalization schemes and plays a role in the 
proof of scheme independence of physical quantities as demonstrated in ref. [jl8 . 



2.5 The Matrices Z^^^ and Z^^^ 

We give here the matrices Z^p^ and Z^^^ in the Feynman gauge. The wave-function renor- 
malization for a quark of charge q is given in the Feynman gauge by 



All € [Att)^ e 2 
where the dots denote terms which are of no interest to us here, and 

Of : 



(2.18) 



2N 



(2.19) 



Applying the renormalization (PTTH) to the quark fields in the operators Qi according 



to 



one finds the matrices Z\ and Z\ . They are given by a single matrix D 



m 



-D, z^ = lcFb. 



(2.20) 



(2.21) 
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The non-vanishing elements of D are 



D(l,l) 


= ^(2,2) 


_ 5 
9 ' 






D(3,3) 


= ^(4,4) 


= ^(5,5) 


= Z)(6,6) 


_ 1 

~ 3 


D(3,9) 


= D(4,10) 


= ^(5,7) 


= ^(6,8) 


_ 2 
9 


D(7,7) 


= ^(8,8) 


= Z)(9,9) 


= Z)(10, 10) 


_ 4 
9 


D(7,5) 


= ^(8,6) 


= ^(9,3) 


= D(10,4) 


_ 1 
9 



3 One-Loop Results 

We give here one-loop results in QED. The corresponding QCD expressions can be found 
in section 3 of ref. [|l^. 

3.1 Current-Current Contributions to 7^*^^ 

The non-vanishing contributions of diagrams Di - D3 of fig. 2 to the matrix 7^°-* are given 
as follows 





cc 


7(0) (2, 2) 


cc 


7f (3,9); 


cc 




_ _8 

cc 3 


7f (5,7); 


cc 


7(°)(6,8y 


cc 


8 

3 ' 








7^(7,5); 


cc 


7i°)(r,7) 


cc 




cc 


;7(°)(8,8); 


_ 4 
3 ' 

cc 


7^(9,3); 


cc 


7(°)(9,9) 


cc 


7(°)(10,4) 


cc 


7^(10,10) 


_ _4 

- cc 



(3.1) 



These results already include the contributions from the wave function renormalization. 



3.2 Penguin Contributions to 7(^) 

The contributions of diagrams Pg^"* ^.nd Pq^^ of fig. 3 to the matrix 7^'^^ have a very simple 
structure. The insertion of any operator of ( |2.1|) into the diagrams of fig. 3 results always 
into the sum of the operators Q7 and Qg multiplied by an overall factor. Denoting by 

P = (0,0,0,0,0,0,1,0,1,0) (3.2) 
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the row vector in the space (Qi, • • • , Qio), the elements of 7^°-* coming from the diagrams 
of fig. 3 are as follows 



7f(Q9) 



27 ' 



^)p 

2N J • 



10) 



16 p 
27 ' 

§{u-l)P. 



N]P. 



B{u + ^]P 



27 



i(. + | + f)P. 



(3.3) 



where u and d {u + d = f) denote the number of effective up- and down-quark flavours, 
respectively. 

It should be stressed that these results do not depend on renormalization scheme for 
Qi and are also valid for Qi. 



3.3 General Comments on f 

As already remarked in section 2.4, the finite terms fg and fg depend on the renormaliza- 
tion scheme used. In the case of the NDR scheme considered here, the penguin diagram 
contributions to these finite terms may also depend on whether Qi or their Fierz conju- 
gates Qi are inserted in the penguin diagrams. This observation has already been made in 
ref. Ijl^ and has been used to avoid a direct evaluation of the two-loop diagrams involving 



Tr (75 7^ 7j,7^7o-), i.e. the type 1 penguin insertions. Here we only state that whereas the 
results for (V — A) ® (V + A) operators Q^-Qs do not depend on the form used, this is 
no longer the case for the {V — A) ® {V — A) operators Qi-Qi, Qq, and Qiq. We will 
return to this in section 5.1, where a generalization of the method of ref. to the mixed 
QED-QCD case is presented. 

4 Current-Current Contributions to the 0{a a^) Ano- 
malous Dimension Matrix j^^^ 

The calculation of current-current contributions to the 0{aas) anomalous dimension ma- 
trix 7£) uses the singularities found in our previous calculation of the C(«s) contributions. 
Due to the fact that certain diagrams present in the pure QCD case do not contribute here 
and due to different colour and electric charge factors, the basic structure of the 0{aas) 
results is quite different from pure QCD. In particular, the mixing between different opera- 
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tors can be divided into three blocks {Qi, Q2), {Q3, Qa, Q9, Qio) and {Q5, Qe, Q7, Qs) with 
no mixing between different blocks. 

The results in this section include the contribution of the matrix Z^P . 



4.1 {V -A)^{V - A) Operators 

The 2x2 matrix describing the mixing in the current-current sector {Qi, Q2) is given by 



8A^ 



22 j_ 

3 N 



25 
3 



_ 2 
3 

3 N 



(4.1) 



The 4x4 matrix describing the mixing in the current-current sector (Qs, Qi,Q9,Qio) 
is given by 

/ 4 -6 RM-^l ^ \ 



6_ 

N 

-3 

AM - 20 J. 
^■^^ 3 N 

11 
3 



3 N 



3 N 
34 
3 

^-'^ 3 N 



16 
3 

^-'^ 3 N 
_ 10 
3 

_2A^- 2i 



(4.2) 



4.2 (y - A) (g) (y + ^) Operators 

The 4x4 matrix describing the mixing in the current-current sector {Q^, Qq, Qy, Qg) is 
given by 

/ -f 6 -m-'^l 4? \ 



^-'^ 3iV 



v 



3N-§ 

20 

3 

3 3 Af 



3 JV 



4/V - 26 J_ 



40 
3 

22 AT _ 16 J_ 

3 3 N 

38 

3 

20 AT _ 26 J_ 
3 3 N 



(4.3) 



3 3 Af 

We observe that the elements of both matrices grow at most as in the Large-A?^ limit 
and consequently the anomalous dimension matrix 0{a ctg) resulting from current-current 
diagrams approaches a constant A^-independent matrix. 



5 Penguin Diagram Contributions to the 0{aas) Ano- 
malous Dimension Matrix 

5.1 General Structure 

The calculation of the penguin diagram contributions to can be considerably simplified 
by first analyzing the general structure of these contributions. The two-loop penguin 
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diagrams are shown in fig. 5 where two types of insertions of a given operator, type 1 and 
type 2, have to be considered. It is known that type 1 insertions are problematic in the NDR 
scheme in which closed fermion loops involving 75 can not be calculated unambiguously. 
Yet, as we have demonstrated in our previous paper, the calculation of 0(a;g) anomalous 
dimensions can be reduced to the calculation of type 2 insertions only. Type 2 insertions do 
not pose any problems. Indeed with the help of eq. (|2.17| ) with a and h standing this time 



for two different bases (see eqs. (|2.1| ) and (|2.2| )) it is possible to relate the type 1 insertions 
to the type 2 insertions. This procedure can be generalized to 0{aas) so that also in this 
case explicit calculations of closed fermion loops can be avoided. The procedure this time is 
more complicated because the insertions of operators involving uu have to be distinguished 
from the insertions of operators involving dd. We will now explain the procedure in some 
detail. 

In order to solve the problem of closed fermion loops involving 75 in the NDR scheme, 
we have to consider four different bases of operators. All four bases contain the penguin 
operators Q3 - Qiq of but differ in the current-current operators Qi and Q2- The 
latter are given as follows. 



Basis A: 



Basis B: 



Basis C: 



Basis D: 



g^"^ = {S^Uli)y_j^{Ul3d^)y_j^ , Q^2^ = {SU)y^^{ud)y^^ . (S.l) 



Q? = (M/3)y_^ (rf/3rfa)^^_^ , = {sd)y _j^{dd) ^ _^ . (5.2) 



Qi") = {sd)y_^{uu)y_j^ , Q^2^ = {So,dp)y_^{upUc,)y_^ . (5.3) 



Qi^ = i.sd)y_j^{dd)^_^, Q^^^ = {so,dp)y_^(dpd^)^_^. (5.4) 

The basis A is the standard basis of ( p.l|) and the basis B is an auxiliary basis needed 
for the solution of the problem. The bases C and D are simply Fierz conjugates of Qi and 
Q2 in A and B, respectively. Evidently, 

Qf^ = Qf^ , Q?^ = Qi'^ . (5.5) 

Let us next denote by [Qi\i and [Qi\2 the result of the type 1 and type 2 insertions of 
an operator Qj, respectively. Then the results for penguin contributions to the row entries 



10 







1 + d [Qf^] 


1 + 2 [Qf^h 






1 + d [q(^)] 


1 + 2 [gf ]2 


[Q9I 


= u[Qt^ 




]l-[Qi'^]2, 


[Qio]p 




1 - [Q? 


]l-[Qf]2. 



of 7g^g) for Q3, Q4, Qg, and Qio can be written as follows 

(5.6) 
(5.7) 
(5.8) 

(5.9) 

Evidently, in order to calculate these rows, the problem of closed fermion loops has to 
be considered. [Qi]p and [Q2]p receive contributions only from type 2 insertions, i.e. 

[Ql]p = [Qi"^]2, [Q2]p = [Q^2^]2, (5.10) 

and do not pose any problems. The case of {V — A) ® {V + A) operators will be discussed 
below. 

There are six independent entries in eqs. ( |5.6| ) - ( |5.9| ), which have to be found in order 
to complete the calculation. We show how this can be done in four steps: 
Step 1: 

[Q[% and [Qi% can be calculated without any problems. The results are given in 
section 5.2. Moreover one finds the relation 

02^ = -^[Qh2. (5.11) 

Step 2: 

Q^^^ and Q^2^ receive contributions only from type 1 insertions. This allows to find 
[Qi"^]i ciiid [<52"'']i by comparing the two-loop anomalous dimension matrices calculated in 
the bases C and A using the relation 



m',' = a^e^r + [(^^s)p,7f J + [(Afe)„7f^J , (5.12) 

where this time 

Because the insertions in the current-current diagrams are identical for these two bases, 
only penguin diagram contributions to j^l^ and Af j, i = s,e enter this relation. On the 
other hand 7^°-* and 'j^'^^ are full one-loop matrices. A simple calculation of finite terms in 
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one-loop penguin diagrams gives 



(Afe)p 



8 

27 



P 




/ \ 

P 


V ) 



(5.14) 



with P defined in (|3.2| ) and 

P -- 



0, 0, 



1 

iV' 



1, - 



1, 0, 0, 0, 



(5.15) 



as in ref. [|T^. The evaluation of ( |5.14| ) does not involve the dangerous traces with 75. 
Step 3: 

[Qi''^]! can be related to by inspecting the diagrams of fig. |. We find 



(5.16) 



and consequently can be found by using ( p.l6| ) and [Qi""*]! obtained in step 2. 

Step 4: 

The calculation of [(52'^''] 1 is slightly more complicated because the inspection of the 
diagrams of fig. ^ does not allow for a simple relation like ( |5.16| ). Since Q^f* receives 
contributions from both type 1 and 2 insertions we can write 



&2\ = [Qi% 



[Q 



2 J2 , 



(5.17) 



with the last entry calculated in step 1. 

In order to find [(52'^^]p we compare the two-loop anomalous dimension matrices calcu- 
lated in the bases D and B using the relation 



(5.18) 



where the index d indicates that now the auxiliary bases D and B are considered. We find 
first 



(Are)i'^) 



— (iV- 11 
27 ^ ' 



( P \ 

-P 




(Ars)S,^) 



/ P \ 

-P 




(5.19) 
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The matrix {'j^^^)^^^ differs from 'j^^'' only in the first row which in this case equals the 
second row. An explicit formula for 7^^'') can be found in appendix A of ref. [15|. The 
matrix (7^*^^ ) is given by eqs. 



, except for the following changes, 



id) 



7(0) (2, 2) 



(d) 



4 
3 



:,(0)(gW. 



27 



iN + l)P. 



(5.20) 
(5.21) 



= 7e ' IV2 

In this way, the commutators in ( |5.18|) can be calculated. Since moreover both types 
of insertions of Qi''^ have been calculated in steps 1 and 3, the element [(52'^^] 1 can finally 
be extracted from ( [5.18|) when in addition the relation ( |5T5| ) is used. 

In summary, we have demonstrated that the penguin contributions to the rows for Qs, 
Q4, Qg, and Qio can be obtained by considering only type 2 insertions. 

The calculation of the rows for Q5 - in the anomalous dimension matrix is simplified 
by the fact that the Fierz symmetry is preserved for these operators in the NDR scheme 
as already discussed in section 3 and ref. [|T^. Thus it is sufficient to consider only the 
operators Q5 and Qq which receive only contributions from type 2 insertions. This gives 
directly the Q5 and Qq lines. In order to obtain the Q7 and Qs lines one has to replace d 
by —d/2 in the Q5 and Qq lines, respectively. 



5.2 Results 



The singular terms in the diagrams of fig. 5 with type 2 insertions can be found in tables 
2-4 of ref. Using these tables, including properly colour and electric charge factors, 
we can find all the insertions necessary to calculate the full matrix (7s^e'*)p according to the 
procedure of the previous subsection. 



We find 



7se,NDR 



243 
+ 

+ 



[N{Q, + Qe) - (Qs + Q5)] 

Qs 
Qio ' 





64- 


Qr 










+ 


I'M 


.9 


27. 


27 




80- 


Q9 










+ 




.9 


27. 


.27 . 



(5.22) 



7ie,NDR(<52) 



556 
243 

+ 



+ 



(Q4 + Qe) 
200 



1 

N 



{Qs + Q5 



243 
200 

243^ 



N + 
N- 



1316 1 ■ 


Qj- 


■124- 


243 A^. 


[ 27 J 


1348 1 ■ 


Q9 + 


■172- 


243 A^. 


[ 27 J 



Qio 



(5.23) 



13 



{u - d/2) — [N{Q, + Q,) - {Qs + Q,)] 



+ {u - d/2) 
+ {u - d/2) 



9 



112 
~27 



-N^ 

9 27 



Q7 + 
Q9^ 



[u - d/2) -N 
[u - d/2) AiV 



(5.24) 



7se,NDR(Q6) 



and 



7se,NDRWl 



{Qa + 



1 



+ 
+ 
+ 



, 64^^ 56 1\ ^/136,, 260 1 
[u - d/2) -—N - — - + / —N + 



27 9 iV 

232 44 1 

f(.-./2,-|/]0, 

,„-./2)(-Pa. ™^ 



V243 



27 ^ ' ' 27 . 



27 iV 
Qio 



+ / 



136 
243' 



243 iV 



100 1 

243 iV 



+ 



+ 



116 

243 
232 



(Q4 + g6)-^(g3 + g5) 



243 243 N 

232 ,^ 920 1 

N A 

243 243 



Q9 



32 

27 
128 



27 



Qio 



Tse, NDr(Q2 ') 



Tse, ndr(Qi) 



(5.25) 



(5.26) 



(5.27) 



where we have indicated that these resuhs have been obtained in the NDR scheme. We 



observe that 



7, 



(0) 



(see eq. (|3.3| )) and ■y^j^ contain elements growing hke N and A^^ 
respectively. Consequently, the corresponding elements in the one-loop 0{a) and two- 
loop 0{aas) matrices grow like N in the Large-iV limit. 



6 Full two— loop Anomalous Dimension Matrix 7, 



(1) 



6.1 Basic Result of this Paper 

Adding the current-current and penguin contributions to 'j^i^ found in sections 4 and 5, 
respectively, we obtain the complete 0{aas) anomalous dimension matrix 7ie,NDR- 

We first note that 16 elements of this matrix vanish. These are the entries [7£NDR]«i 
with i = 3, . . .,10; j = 1,2. 
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The remaining 84 elements are non-vanishing. These entries of 7^0 ndr given in 
table |I]. For phenomenological applications, we need only the results with = 3. We give 
them in appendix C for an arbitrary number of flavours. 

Let us just comment on certain features of the matrix 7fe,NDR- 

• A comparison of 7^°-' and 7fe ndr shows that the impact of QCD is to fill out most 



of the zero entries present in 7 



(0)_ 



The corrections to non-vanishing elements in 7^°^ are moderate. For = 3 and 
u = d = 2 the largest corrections are found in elemets (7,5), (8,6) and (4,7). For 
as/4:TT ^ 0.01 they amount to 10 — 15 % corrections. Thus the corrections are 
certainly smaller than in the case of 7sndr evaluated in ref. [|1^. Similar comments 
apply to elements which become non-zero at 0{aas)- They are smaller than the 
corresponding terms in the matrix 7^^^^ 

It is needless to say that all these comments are specific to the NDR scheme and the 
true size of the next-to-leading order corrections can only be assessed after a full 
renormalization group analysis has been performed. We will return to this question 



in ref. [18 



As in the case of the Ola"^) corrections, on general grounds the coefficient of the 1/e^- 
divergences at 0{aas) in the unrenormalized Green function, f I2, is entirely given in terms 
of quantities calculable at one-loop. 



22 



if if + ifif 



ai7 



(0) 



(6.1) 



with ai and 'j^^^ given in ref. |jT5| . The explicit calculation shows that this relation is indeed 



satisfied, which constitutes a check of our calculation. For details, the reader is referred to 



sect. 6.3 of ref. [15 . 



6.2 Comments on the HV Scheme 



Having the two-loop anomalous dimension matrix 7ie ndr hand, it is a simple matter 
to obtain this matrix in any other scheme by using the relation ( p.l7| ). As an example, we 
consider here the HV scheme for which 7^ hv has been explicitly given in ref. [|T5[ . 
In order to find 7se hv simplest method is to use 



7se,HV ~ /se,NDR ^ 



Afs,7, 



(0) 



+ 



Are, 7 



(0) 



(6.2) 
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where this time 



Afs = (fs)HV - (^s)nDR , Afe = (fe)HV " (^e)NDR • (6.3) 



The result for Afs can be found in sections 3.3 and 3.4 of ref. [0. Here we give 
in addition the result for Afg. As in the case of Afg it is convenient to separate the 
contributions to Afe into current-current and penguin contributions obtained from finite 
terms in figs. |^ and ||, respectively. 

The non-vanishing elements in [Afej^c are as follows: 



[Are(l,l)],e = [Afe(2,2)]^^ = -[Afe(9,3)],, = -[Afe(10,4)]^ 

[Afe(3,3)],, = [Afe(4,4)]^^ = [Afe(5,5)],, = [Afe(6,6)],, 

[Afe(3,9)],, = [Afe(4, 10)],, = [Afe(7,7)],, = [Afe(8,8)],, 

[Afe(9,9)],, = [Afe(10,10)],, = |, 

[Afe(5,7)],, = [Afe(6,8)],, = 2[Afe(7,5)],, = 2 [Afe(8, 6)],, 



For the penguin contributions one gets: 

[Afe(Ql)], = -[Afe(Q4)]p = 2[Afe(gio)]p = -^.N?.^ 



[Afe(g2)]p = -[Afe(g3)]p = 2 [Afe(g9)]p = "^P , 
[Afe(Q5)]p = [Afe(g6)]p = [Afe(Q7)]p = [^UQs\ = 



2 

9 ' 
2 

3 ' 
4 
'9 ' 



20 
■ 9 



(6.4) 



(6.5) 



With the help of eqs. (|6.2| )- (|6.5| ) and Afg of ref. ||T5|, one can find an explicit expression 
for 7£hv- Sive it here since due to the complexity of this matrix it is anyway 

better to make the transformation (|6.2|) by computer. 



7 Summary 

We have presented the details and the explicit results of the calculation of the 10 x 10 
two-loop anomalous dimension matrix 0{aas) involving current-current, QCD-penguin 
and electroweak penguin operators. 

Performing the calculation in the simplest renormalization scheme with anticommuting 
75 (NDR) we have demonstrated that a direct evaluation of penguin diagrams with 75 
appearing in closed fermion loops can be avoided by studying simultaneously four different 
bases of operators. In this way an unambiguous result for 7^^ in the NDR scheme could 
be obtained. 
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This analysis completes our extensive calculations of the anomalous dimension matrix 
7 in eq. (|1 . 1| ) which we have presented in |13, 14, ^ and in the present paper. This matrix 
constitutes an important ingredient of any analysis of non-leptonic weak decays which goes 
beyond the leading logarithmic approximation. The full next -to-leading renormalization 
group analysis of the Wilson coefficient functions of operators Qi has been just completed 
and the details of this work together with phenomenological implications can be found in 
ref. [ITHl. 
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Appendices 

A One— Loop Anomalous Dimension Matrix 7^°^ 



8 
3 

















16 AT 
27 





IGN 
27 








8 
3 














16 
27 





16 
27 























16 1 16 N {u-d/2 ) 
27 27 





88 1 16 Af (u-d/2) 
27 27 























-16 Af 1 16(«-d/2 ) 





-16 AT 1 16 (u-d/2) 


8 


27 ' 27 


27 ' 27 


3 




















8 1 16 AT (u-d/2 ) 
3 27 





16 AT (u-d/2) 
27 























16 {u-d/2) 
27 


8 
3 


16 (u-d/2) 
27 

















4 
3 





4 1 16N{u+d/4 ) 
3 27 





16 AT (u+d/4) 
27 




















4 

3 


16 {u+d/4) 
27 


4 
3 


16 (u+d/4) 
27 











4 
3 











8 1 l6N{u+d/4) 
27 ~'~ 27 





28 1 16 AT (u+d/4) 
27 27 














4 
3 








8 Af 1 16 {u+d/4) 
27 27 





8 Af 1 16 (u+d/4) 
27 27 


4 
3 
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B Table for Two— Loop Anomalous Dimension Matrix 
7^) in the NDR Scheme 



Table 1: Pull QCD-QED Anomalous Dimension Matrix {■yil^)ij for the NDR scheme (with 
vanishing entries omitted tacitly). 



'■J) 



1,1 
1,3 
1,5 
1,7 
1,9 



-f + 87V 



243 
88 
243 

_64 I SAT^ 
27 9 

_80 _|_ 8Ar2 
27 9 



1,2) 
1,4) 
1,6) 
1,8) 
1,10) 



_2 

3 

88 

243 
88 AT 

243 
40 Af 

27 
56 Af 

27 



2,1 
2,3 
2,5 
2,7 
2,9 



25 
3 

-556 
243 AT 
-556 
243 Af 
1316 200 Af 
243 
200 AT 



243 AT 
-1348 
243 Af 



243 



2,2) 
2,4) 
2,6) 
2,8) 
(2, 10) 



3N 

556 
243 
556 
243 
124 
27 
172 
27 



3,3 

3,5 
3,7 



3,9 



-1040 
243 Af 
-1400 
243 Af 



1690 _ 136 {u-d/2) 
243 N 243 

232 _ 136 (u-d/2) 
243 TV 243 
464 AT _ 112 (m -d/2) , 

243 27 
2408 AT _ 32 (m -d/2) _^ 
243 27 



SAT^ {u-d/2) 
9 

8N^ {u-d/2) 



3,4) 
3,6) 
3,8) 
(3, 10) 



243 



243 



1690 _|_ 136 AT {u-d/2) 
' 243 243 

232 I 136 A^ {u-d/2) 

243 243 

64 _|_ 88 N {u-d/2) 

27 ~^ 27 
112 I 8 N {u-d/2) 
27 27 



4,3 
4,5 
4,7 
4.9 



641 I 
243 



243 



+ 



AT 

-388 m i 32 d 
243 "'' 243 



Af 



64 
27 

386 _ 
27 



- 8Af^ J- Af f 280 m I 64d 
9 V 243 ~'~ 243 

8N^ _i_ ]\f ( 280 u I 64_d\ , 
9 V 243 ""^ 243 / 



N 

-1612 11 I 



N 



4,4) 
4,6) 
4,8) 
(4, 10) 



_6_ _ 817 Af _|_ 388m _ 32d 
Af 243 ~'~ 243 243 



-88 Af 1 
243 ~'~ 


388 « 


32 d 


243 


243 


-40 Af 


100 m 


16 d 


27 


27 


27 



-40 _ 2Af _|_ 148m _ 80d 
3Af 27 ~'~ 27 27 



5,3 
5,5 
5,7 
5,9 



-16 

3Af 



-136 {u-d/2) 
243 

-6 _ 136 (M-d/2) 
Af 243 
g ^ _ 112 {u-d/2) _^ 8Af^ {u-d/2) 

J2 , 



27 



-32 {u-d/2) I 8N^{u -d/2) 
27 9 



5,4) 
5,6) 
5,8) 
(5, 10) 



136Af (M-d/2) 
243 

^ I 136 Af (M-d/2) 
" 243 
40 I 88Af (M-d/2) 
3 27 
8 N {u-d/2) 
27 
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Table 1: Full QCD-QED Anomalous Dimension Matrix for the NDR scheme (con- 

tinued; with vanishing entries omitted tacitly). 



[hJ) 



■21-2 d 
243 



(6,3) 
(6,5) 



N 

-748 u 



-2 + AT ( 



3 + 

-440 u 



212 d 
243 



243 

1000 d 



+ 



N 



+ 



N 
424 d 
243 



-440 u 
243 



+ 



N 
424 d \ 
243 / 



(6,4 
(6,6 
(6,8 
(6,10) 



-6 
N 
-16 
3N 



+ 



748 u 

243 



22 N 



212 d 




243 




748 m 


212 d 


243 


243 


1 188 m 


160 d 


' 27 


27 



-140m _|_ 64d 
27 "'~ 27 



-136 {u+d/4) 
243 



— - 4 

3iV ^-^^ 



136 {u+d/4) 
24c 

+ 



3 AT 243 

112 (M-l-rf/4) , SAT^ {u+d/4) 
27 ~'~ 9 



-32 (M+(i/4) , 8N^ {u +d/4) 

27 9 



(7,4 
(7,6 
(7,8 

(7, 10) 



136 N {u+d/4) 
243 



+ 



136 N {u+d/4) 
243 



20 
3 

38 I 88 AT {u+d/4) 
3 



SAf (m- 



27 

-d/4) 



243 



27 



.3) 
5) 
7) 

,9) 



243 



2 + 



243 



A^ 



N 
-440 m 
243 



-1 + 

508 d 
243 



N 

-748 u 
243 



106 d 
243 



+ 

212 d 

243 



N 
+ 



-440 m 
243 

170011 
243 



N 



_ 212 d 

243 

500 d 
243 



(8,4 
(8,6 
(8,8 
(8,10) 



748 m _|_ 106 



243 



243 



-8 I 11 Af _|_ 748 m _|_ 

3Af 3 243 ~'~ 

-26 I 20 Af I 188 m , 

3Ar 3 27 



106 d 
243 
80 d 
27 



-140 m 
27 



32 d 
27 



243 AT ^ 



520 
243 A/" ■ 
-1082 
243 Af 



+ 



-116 
243 A^ 
232 AT 
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740 AT 



136 (M+d/4) 
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243 



112 (m+ d/4) I 8 (M-l-d/4) 
27 9 



32 (m 



243 



d/4) , 8N^ (M-l-d/4) 
27 9 



(9,4 
(9,6 
(9,8 
(9,10) 



243 
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243 



764 _|_ 
243 ~'~ 
116 _|_ 
243 ~'~ 
_32 I 
27 

38 I 
27 



136 N (m4 


d/4) 


243 




136 N (m4 


d/4) 


243 




88Af (m4 


-d/4) 



SAr(M-|-d/4) 
27 



:io,3) 
:io,5) 

;io,7) 
;io,9) 



32 
' 27 



32 I 
27 9 



I 4N^ 
^ 9 
4Af2 



+ 



1333 

243 
44 
243 

620 m 
243 



+ ■ 
+ 



AT 
-388 m 
243 



40 d 
" 243 



N 



16 d 
243 



+ 



N 

-1612 u 328 d 
243 



243 



N 



+ N ( 

+ N ( 



280 M 
243 

280 M 
243 



32 d 
243 
32 d 
243 



(10,4) 
(10,6) 

(10,8) 
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-20 _|_ 287 Af I 388m , 16d 
3N ^ 243 243 243 

44 Af I 388 m I 16 d 
243 ~'~ 243 ~'~ 243 



100 u 



27 



20 Af _ 
27 

-2 _ 26 Af 

3Af 27 ~'~ 27 



8d 
27 



148 m I 40d 
27 
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C Two— Loop QCD— QED Anomalous Dimension Ma- 
trix 7g^g) for = 3 in the NDR Scheme 
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212 d 
243 



136(M+d/4 ) 
81 
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243 

556 
729 

232 _ 136 {u-d/2) 
729 243 



243 



388m _|_ 32d 
729 729 



r, _ 136 (M-d/2) 
243 

748 m _|_ 212 d 



729 



729 



116 
9 

-1 - 

116 
729 

44 

243 " 



136 (M+d/4) 
243 



748 m 
729 



106 d 
729 



136 (M+d/4) 
243 



388 m 
729 



16 d 
729 



88 152 40 136 56 

81 27 9 27 9 

556 _484 _124 3148 172 

243 729 27 729 27 

232 I 136 (M-d/2) 3136 , 104 (M-d/2) 64 , 88 (M-d/2) 20272 , 184 (M-d/2) _112 , 8 (M-d/2) 

243 ~'~ 81 729 ^ 27 27 ~'~ 9 729 27 27 9 

_88 _|_ 388 m _ 32d _152 , 3140 m , 656 d _40 _ 100 m _ 16d 170 , 908 m , 1232 d _14 , 148 m _ 80d 

81 243 243 27 729 729 9 27 27 27 729 729 3 27 27 

r. I 136 (M-d/2) _232 I 104 {u-d/2) 40 , 88 (m -d/2) 184 (M-d/2) 8 (M-d/2) 

" 81 9 27 3 9 27 9 

7 _l_ 748m _ 212 d o _ 5212m , 4832 d 182 , 188 m _ 160 d -2260m , 2816 d -140m , 64d 
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20 _|_ 136 (M+d/4) 134 , 104 (M+d/4) 38 , 88 ( M+d/4) 184 (M+d/4) 8 (M+d/4) 

3 81 9 27 3 9 27 9 
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Figures of Feynman Diagrams 




(1) (2) 

Figure 2: 




(1) (2) 



Figure 3: 
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Figure 4: 
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(1) 



(2) 



(3) 



(4) 




Figure 5: 
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Figure Captions 



Figure 1: The three basic ways of inserting a given operator into a four-point function: 
(a) current-current-, (b) type 1 penguin-, (c) type 2 penguin-insertion. The wavy 
hues denote gluons or photons. The 4-vertices "® ®" denote standard operator 
insertions. 

Figure 2: One-loop current-current diagrams contributing to 7^°\ The meaning of hues 
and vertices is the same as in fig. 1. Possible left-right or up-down reflected diagrams 
are not shown. 

Figure 3: One-loop type 1 and 2 penguin diagrams contributing to 7^°-'. The meaning of 
lines and vertices is the same as in fig. 1. 

Figure 4: Two-loop current-current diagrams contributing to The meaning of lines 
and vertices is the same as in fig. 1. Possible left-right or up-down reflected diagrams 
are not shown. 

Figure 5: Two-loop penguin diagrams contributing to The wavy lines denote gluons 
or photons. Square- vert ices stand for type 1 and 2 penguin insertions as of figs. ||(b) 
and (c), respectively. Possible left-right reflected diagrams are not shown. 
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